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ABSTRACT 

There is a robust upper limit on the energy of synchrotron radiation in high-energy astro- 
physics: ~ moC^/a, where a = 1/137 is the fine structure constant and the value refers to the 
comoving frame of the fluid. This is the maximal energy of synchrotron photons which can 
be emitted by an electron having an arbitrarily high initial energy after it turns by angle ^ tt 
in the magnetic field. This upper limit can be naturally reached if the converter mechanism 
contributes to the jet radiation and can be imprinted in spectra of some blazars as a cutoff or a 
dip in the GeV range. We use numerical simulations to probe the range of parameters of a ra- 
diating jet where the ultimate synchrotron cutoff appears. We reproduce the variety of spectra 
depending on the source luminosity and on the scale of the emission site. We also compare our 
results with the EGRET blazar spectra in order to illustrate that agreement is possible but still 
not statistically significant. The predicted feature, if it exists, should be observed by Fermi in 
spectra of some blazars. 



1 INTRODUCTION 

The issue of maximal possible synchrotron radiation energy in as- 
, trophysics (Guilbert, Fabian & Rees 1983) was discussed in var- 
ious frameworks: the extreme Fermi acceleration of electrons at 
shocks (e.g. Ghisellini 1999) and the converter mechanism (Der- 
ishev 2003; Stern 2003). In the extreme acceleration scenario the 
' electron energy is limited by the condition that the Larmor gyra- 
, tion time, being the shortest available acceleration time, is shorter 
than the synchrotron cooling time. In the converter mechanism, a 
new high-energy electron appears in the relativistic jet as a result 
of photon-photon pair (or photomeson) production and its initial 
energy can be arbitrarily high in the fluid frame. However, its ini- 
tial direction in the jet frame is opposite to the jet bulk motion. 
Therefore, its synchrotron radiation at the first moment is Doppler 
deboosted in the external reference frame. When the electron has 
turned around, the emitted synchrotron photons become Doppler 
boosted. When it gyrates in the magnetic field by the angle vr, 
it cools down below the limiting Lorentz factor (in the comov- 
ing frame) 7m ~ 10*5"^''^, where B[G\ is the magnetic field 
strength. In the extreme acceleration scenario it is approximately 
the same and in both cases the maximal comoving energy of syn- 
chrotron photons does not depend on the magnetic field: Sm ~ 1/q 
(i.e. ~ lO^MeV), hereafter e — E/rricC? is the photon energy 
in electron rest mass units. In the case of blazars, this cutoff is 
blueshifted by factor up to 2r, where F is the bulk Lorentz fac- 
tor of the jet. Observations of the superluminal proper motion in 
blazars imply that a typical Lorentz factor is about F ~ 10 — 20 
(see Jorstad et al. 2001). In some cases F ~ 40 and its value in 
the emission region could be even higher, when taking into account 
possible deceleration of the jet (see Stem & Poutanen 2008, here- 
after SP08). Therefore, the ultimate synchrotron cutoff in the blazar 
spectra should be expected in the GeV range, maybe slightly above 



10 GeV in some extreme cases. This feature should be seen as a 
spectral dip with a continuation of the spectrum to higher energies 
as a result of Comptonization of the external or synchrotron radia- 
tion. 

In which physical context can the observable synchrotron cut- 
off in GeV range appear? The main mechanism of particle acceler- 
ation, which is responsible for the gamma-ray emission of blazars 
is still unknown. There are at least three possibilities which are dis- 
cussed: 

(i) Diffusive acceleration of electrons by various plasma per- 
turbations or turbulence (Fermi II mechanism). The ultimate ac- 
celeration in this case would imply that the relative bulk veloc- 
ities of the fluid at the gyroradius scale are relativistic (in the 
conditions of examples considered here the maximal gyroradius 
Tui ~ 10~^ — 10~'*i?j, where _Rj is the jet radius). This is scarcely 
possible. 

(ii) Shock acceleration of charged particles (Fermi I mecha- 
nism). It can work at internal shocks in the jet or at the shear layer 
between the jet and the external environment. In this case, acceler- 
ation rate can be much faster, especially at the shear layer, where 
an electron entering the jet from the external environment can gain 
factor in energy (provided that the thickness of the transition 
layer is less than the gyroradius). Therefore, the particle spectra 
can extend up to 7m- The issue is how much power can be emitted 
by the particles near the cutoff. Simulations of shock acceleration 
(Niemiec & Ostrowski 2006; Ostrowski 2008) demonstrate a large 
variety of resulting particle spectra depending on the geometry of 
the magnetic field and the power spectrum of its perturbations. For 
some parameters the spectra are very hard. On the other hand, in 
these simulations the radiative energy losses are neglected and it is 
difficult to conclude whether they extend close to 7m- In the shear 
layer scenario, an electron having the energy above 7m cannot pen- 
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etrate into the jet deeper than rm- Therefore, the energy budget for 
extreme acceleration at the shear layer is restricted by the bulk en- 
ergy of a thin outer layer of the jet. 

(iii) Converter mechanism, which can take a form of a run- 
away photon breeding in relativistic jets (Derishev et al. 2003; Stern 
2003; Stern & Poutanen 2006). The photon breeding mechanism 
is probably the most efficient way of dissipation of the relativis- 
tic jet energy into gamma-rays. The mechanism should work if 
the active galactic nucleus (AGN) is sufficiently powerful (above 
Ld ~ 10''^ — lO^'^erg s~^, depending on the scale of the emis- 
sion site) and the if jet is sufficiently fast (the Lorentz factor above 
r ~ 10, see SP08). When photon breeding mechanism operates, 
the instant energy gain at photon - electron-positron conversion 
in the jet can reach AT^ and the ultimate energy 7m can be eas- 
ily reached. Stem & Poutanen (2006), studying the photon breed- 
ing regime in jets, showed that the energy of particles through 
the breeding cycle grows till the Lorentz factor of produced pairs 
reaches 7m in the jet frame. SP08 performing numerical simula- 
tions of the photon breeding have obtained some hard (/3 ~ 2) 
spectra at a moderate luminosity of the source, where the ultimate 
synchrotron cutoff is clearly seen (hereafter we define the photon 
spectral index (3 as dN/de oc £~^). As discussed above, the possi- 
bility to reproduce such cutoff with other mechanisms is uncertain 
at best. Therefore, the observation of such a feature would be an 
argument in favor of the converter mechanism as a source of the jet 
radiation. 

Existing high-energy data supplied mainly by EGRET (30 
MeV - 20 GeV energy range), do not reveal such a cutoff at a 
significant level. Nevertheless, there are some indications of a non- 
trivial shape of the blazar spectra in the high-energy range, more 
complicated than a single high-energy hump. Such an indication 
has been found by Costamante, Aharonian & Khangulyan (2007) in 
the spectrum of BL Lac PKS 2155-304. This is one of five objects, 
which has been detected by both EGRET and Cherenkov arrays: 
the combination of these two data sets imply that there is a double 
humped structure in the ten MeV - TeV range with a depression 
around 10 GeV. 

This study is motivated by the launch of Fermi (former 
GLAST). This instrument covering 200 MeV - 200 GeV range is 
perfect for detection of the feature we study (see also discussion in 
Costamante et al. 2007). 

In this work we perform a series of simulations in order to out- 
line the variety of conditions, where the ultimate synchrotron cutoff 
can be observed, and to demonstrate the variety of spectral shapes 
in the MeV - TeV range produced by the photon breeding mecha- 
nism. The problem setup, described in Section 2, is similar to that 
of SP08. In the same section we present results of numerical simu- 
lations and discuss them. In Section 3 we compare our results with 
the EGRET blazar data in order to illustrate that the existing data 
have slightly insufficient accuracy to make significant conclusions 
about the structure of blazar spectra in the GeV range. In section 4 
we discuss the perspectives of observations. 



2 THE NUMERICAL SIMULATIONS 

We have performed Monte-Carlo simulations of the jet radiation 
in the photon breeding regime using the Large Particle Monte- 
Carlo Code (LPMC) for particle propagation and interactions and a 
two dimensional ballistic model for the jet dynamics. The problem 
setup and the technique of numerical simulations are described in 



SP08, with details given by Stem & Poutanen (2006). LPMC code 
is described in Stern et al. (1995). 

Briefly, the problem setup can be summarized as follows: 

(i) The jet is represented as a cylinder of radius iij and length 
= 20_Rj split into array 64x20 of concentric shells and slices, 

where each cell can independently decelerate by exchanging 4- 
momentum with interacting photons. The intemal pressure is ne- 
glected. The initial Lorentz factor V at the inlet is constant. The 
magnetic field has toroidal geometry. 

(ii) Initial state. A simulation starts with a uniform jet of a con- 
stant Lorentz factor and a small amount of seed high-energy pho- 
tons. 

(iii) The external environment is kept at rest (see SP08 for the 
corresponding condition) and has the magnetic field perpendicular 
to the jet with random orientation in azimuthal direction. 

(iv) Emission sites. SP08 studied three possible emission sites: 
at the scale of the broad emission line region (R ~ lO^^cm, model 
A), the parsec scale (model B) and small scale (ii ~ 30 — lOO-Rg ~ 
10^* cm, where Rg is the gravitational radius, model C). Here we 
consider models A and C. 

(v) The external radiation. For model A we assume two- 
component radiation field, which includes direct radiation of the ac- 
cretion disc (the multicolor blackbody spectmm and the X-ray tail) 
and the isotropic radiation with the spectrum extending from the 
UV (lines and scattered disc radiation) to the far IR (dust radiation 
from a few parsec scale and a possible contribution of synchrotron 
radiation). The energy density of the isotropic component is 0.05 
of that of the direct disc radiation. The spectrum of the isotropic 
component is parametrised as a power law dA*'/d£ cx e~^'* in the 
range 10^* < e < Tmax, where Tmax = 0.5 x 10~^ is the maxi- 
mal temperature of the accretion disc in m^t? units. Model C rep- 
resents a smaller distance scale: the vicinity of the accretion disc 
at tens ~ hundreds gravitational radii. At such a distance the pho- 
ton breeding does not require a separate isotropic component of 
the soft background radiation, because the accretion disc supplies 
some photons at sufficiently large angles. In the case of model C, 
the picture of radiation field is more certain than that in model A: 
the temperature of the disc varies as T = rmax(»'/''min)~'^'''', the 
luminosity is distributed as dL/dr oc r~^, the angular distribu- 
tion of the emitted photons flux is described by the Lambert law 
dF/d cos9 oc cos 6. The disc radiation field was precomputed 
as a function of the distance from the central engine in the range 
(30 - 100)i?g. 

In both models, the spectrum of the disc radiation includes 
a nonthermal X-ray tail described by the power law with photon 
index (3 — 2 extending up to 100 keV with the total luminosity 
0.07 of that of the thermal component. Such X-ray emission exists 
in many AGNs. For model C the X-ray component is important, 
because it provides a sufficient opacity for photons of energy e ~ 
10* which are otherwise too soft to interact with the thermal disc 
photons (see Fig.[TJ. 

The main difference between the two models can be expressed 
in terms of dependence of the photon-photon opacity on the direc- 
tion. Fig.[T]shows the opacity for photons moving along the jet and 
in the opposite direction. One can see that in the case of model C, 
the anisotropy in the opacity is higher than that in model A. 

We have chosen these two models as idealised extreme ver- 
sions of extemal radiation: the emission of a point-like disc plus 
an isotropic background (model A) and the radiation of a geomet- 
rically large accretion disc only (model C). Actually the soft ra- 
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Figure 1. The photon-photon opacity at the photon path Rj depending 
on the energy of photons moving in the counter-jet (upper curves for both 
models) and the jet direction (lower curves for both models). The absolute 
value of the opacity corresponds to the disc luminosity 2 X lO^'^erg 
for model A and 2.5 X 10*"^ erg s" 



for model C. 



Figure 2. Spectra produced by the photon breeding mechanism obtained 
with the numerical simulations. Model A: curves from bottom to top coiTe- 
spond to runs A1-A4 (disc luminosity varies from0.5xl0''''crgs-i to2x 
10"'''orgs-i). Model C: curves from bottom to top correspond to runs Cl- 
C6 (the disc luminosity varies from 10''^ erg to 3.2 X 10'*'' erg s~^). 



diation field should be their superposition. Particularly, model C 
should also include some scattered radiation, which we neglect. 

The relevance of each model for different classes of objects 
is not evident. For the operation of the photon breeding in model 
C, a lower luminosity of the accretion disc is required than in 
model A: Ld ~ W^^Riserg versus ~ 10*^i?i7erg s"^ 
(R„ — i?/10"cm). In units of compactness (which is proportional 
to Ld/R), model A has a lower threshold. Model C is probably the 
only one which can work in the case of BL Lacs. On the other hand, 
model C implies that the acceleration of the jet is fast: the Lorentz 
factor r > 20 should be reached at a few tens Rg. Otherwise this 
model does not provide the conditions for the exponential photon 
breeding process. In model A the conversion of the jet energy into 
the radiation is more efficient. 

Other parameters for model A: F = 20, the average distance 
from the black hole 2 x lO^'^ cm, the jet radius Rj — 10^® cm, the 
jet power Lj = Ld, the magnetic energy flux Lb = 0.2Ld- Note 
that the jet power is not important unless nonlinear effects become 
substantial. The latter take place at Lj ~ 10^^ erg s^^. 

Parameters for model C: the distance between the beginning 
of the active region and the black hole R = 10^^ cm, gravitational 
radius (which defines the accretion disc scale) 7?g = 3 x lO^^cm, 
the jet radius Rj = 10^'^cm, Lb = 0.25Z/d, Lj = Ld. 

We have performed a series of simulations for each model 
varying the density of the external soft photon field. This was 
parametrised through the corresponding luminosity of the accre- 
tion disc Ld. In both cases we have started from the luminosity 
which is slightly above the threshold for the supercritical photon 
breeding and have increased it by steps small enough to trace the 
evolution of the spectrum. For model A the simulations (A1-A4) 
correspond to Ld = (0.5, 1, 1.4, 2) x 10*'* erg s^*. For the model 
C (runs C1-C6): Ld = (1, 2, 4, 8, 16, 32) x 10*^ erg s"*. 

The resulting efficiency of the jet power conversion into radi- 
ation varies from 0.16 (run AI) to 0.6 (run A4) for model A and 
from 0.05 (run CI) to 0.3 (run C6) for model C. The resulting spec- 



tra as a function of the disc luminosity for each model are given in 
FigEl 

In both cases the ultimate synchrotron component is relatively 
strong and the cutoff is well seen at a low disc luminosity: within 
factor ~ 3 above the threshold of supercritical photon breeding in 
the case of model A, and within one order of magnitude in the case 
of model C. In real sources, the compactness parameter is proba- 
bly widely distributed, and therefore the number of objects demon- 
strating the ultimate synchrotron cutoff in its clear form should be 
relatively small. At a higher disc luminosity the cutoff turns into a 
shallow dip. In both cases, the feature is associated with the transi- 
tion between the synchrotron and the Compton components, how- 
ever, the position and the shape of the break are different for the 
two models. 

The reasons why the transition between the synchrotron and 
the Compton components shifts to a lower energy as the compact- 
ness increases are as follows: 

(i) In the case of model C the high-energy photons (e ~ 
l/Lmax produced in the external environment are removed from 
the breeding cycle. Indeed, the photon-photon opacity across the 
jet for e — l/Tmax can be as large as a few hundred and most 
of the external high-energy photons are absorbed beyond the jet. 
Only the photons with e < 5 x 10* can penetrate the jet deeply 
enough (see Fig[T) to support efficiently the breeding cycle. The 
synchrotron component induced by these photons is cut off at 

where Bo = 4.4 x 10*''G. Indeed, one can see a maximum at 
e ~ 10 for model C at a high compactness. 

(ii) For model A the high-energy photons produced in the jet 
are absorbed because of a high opacity produced by the isotropic 
soft photons. One can see in Fig[T]that for model A the opacity 
of the isotropic background at the path length ~ WRj exceeds 
unity when e > lO*". The corresponding absorption turnover is seen 
in Fig. [2] In this case, the absorption initiates the electromagnetic 
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cascade inside the jet producing softer pairs and broader spectra 
without sharp features. 

The variety of spectra we present here is incomplete. For ex- 
ample, when the ratio between the magnetic energy flux and the 
disc luminosity is increased, the Compton peak becomes lower and 
the jet radiative efficiency decreases. A larger series of simulations 
will be presented in a separate publication. We also do not discuss 
the observed low-energy synchrotron maximum treated by SP08 
as a secondary effect: radiation of pairs in heating/cooling balance 
further downstream along the jet. 



3 EGRET DATA 

The EGRET catalog (Hartman et al. 1999) includes 66 high- 
confidence identifications of blazars (BL Lac objects, flat-spectrum 
radio quasars, or unidentified flat-spectrum radio sources). 

The data errors are rather large and all spectra of 66 blazars 
with a few exceptions can be fit with a simple power law at a sat- 
isfactory x^- In order to illustrate the present accuracy of observa- 
tions we show four of seven spectra with largest deviations from 
a power law (reduced > 2) in comparison with our simulated 
spectra applying their measured redshift (Fig.[3]l. In these cases our 
spectra give a better description than a power law. In the case of 
3EG 10210-5055 the cutoff in the simulated curve is at a slightly 
smaller energy than is indicated by the data, but this could be com- 
pensated by a larger Lorentz factor of the jet. 

Of course, none of these examples confirms the existence of 
the synchrotron turnover in the GeV range: the statistical signifi- 
cance is insufficient. 

An additional evidence of the high-energy synchrotron maxi- 
mum could be obtained from comparison of the EGRET and the 
TeV spectra for those objects, where both kinds of data exist. 
At present, there are five blazars detected both by EGRET and 
Cherenkov telescopes. Note that observations in the EGRET and 
the TeV ranges are not simultaneous. 

1. PKS 2155-304. As we already mentioned in the Introduc- 
tion this case has been studied by Costamante et al. (2007) as a pos- 
sible example of the synchrotron maximum in the EGRET range. 

2. Mkn 421. It has a comparatively hard spectrum in the 
EGRET range (/3 ~ 1.6) extending up to 10 GeV, which is in- 
consistent with the high-energy synchrotron component unless the 
Lorentz factor of the jet is very large F > 50. 

3. lES 1959-1-650. The TeV spectrum is pubhshed by Aharo- 
nian et al. (2003). The EGRET spectrum is soft (/3 ~ 2.5) with 
the GeV energy flux e^diV/de ~ 10"" TeV s~^ cm"^ while the 
TeV energy flux corrected for the infrared background absorption 
is slightly higher. This resembles the case of PKS 2155-304 and 
intermediate curves in Fig.2. 

4. 3C 279. In the EGRET range the spectrum has a smooth 
maximum at 1 GeV. The spectrum of TeV emission discovered by 
Teshima et al. (2007) has not been published yet. 

5. BL Lacertae. The spectrum is very variable both in the 
EGRET and the TeV ranges (see Albert et al. 2007), and no clear 
conclusions can be made because of the lack of simultaneous ob- 
servations. 

The above consideration just illustrates the present observa- 
tional status. The data do not reject the hypothesis of the existence 
of the ultimate synchrotron cutoff in some blazar spectra, but can- 
not confirm it either, because of the insufficient gamma-ray statis- 
tics. However, it is clear that even a moderate increase in the statis- 
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Figure 3. Examples of the EGRET spectra sampled from the seven spectra 
with largest reduced of a power-law fit resembling some two-component 
simulated curves from Fig. |2] None of these examples gives a significant 
confirmation of the two-componentness because of the large data errors. 
The simulated curves correspond to: run CI at 2=1.003 (3EG J0210-5055), 
the run C4 at 2=0.902 (3EG 1733-1313), the run A2 at 2=0.894 (3EG 
J0540-4402) and the run A2 at 2=1.494 (3EG J1409-0745). 



tics will provide sufficient significance in cases like those presented 
in Fig.[3l 



4 CONCLUSIONS AND PERSPECTIVES OF 
OBSERVATIONS 

Our simulations demonstrate that the ultimate synchrotron cutoff 
in its clear form is a relatively rare phenomenon. In the framework 
of the photon breeding mechanism, the cutoff in the GeV range ap- 
pears when the compactness of the external radiation is just above 
the threshold of the supercritical photon breeding. When the com- 
pactness exceeds the threshold by an order of magnitude, the spec- 
trum becomes almost flat or the synchrotron component moves to 
the softer range (see Fig. |2ll. Another tendency, which can be de- 
rived from the results, is that the spectral sequence depends on the 
angular distribution of the external radiation. When the large-angle 
component of the external radiation in the emission region is rel- 
atively strong (e.g. ~ 0.1 of the direct disc radiation as in model 
A) the spectra are smoother and the separation between the syn- 
chrotron and Compton components is weakly expressed. When the 
side-on radiation is weak, e.g. it is supplied only by the periphery 
of the accretion disc, all features are much sharper. 

The number of blazars with spectra demonstrating the clear 
ultimate cutoff should be relatively small. However, Fermi will be 
able to reveal this kind of spectra if it exists, because it will ob- 
serve simultaneously a large number of objects. The cutoff in the 
right place is by itself insufficient to conclude that we observe the 
ultimate synchrotron cutoff: the latter can also be caused by the 
photon-photon opacity in the same energy range (see SP08). To be 
sure that the observed cutoff has a synchrotron nature, one have 
to observe the continuation of the spectrum above the cutoff in a 
form of the higher energy Compton component extending to the 
TeV range. In principle, Fermi can detect the low-energy part of 
the Compton maximum at i5 '--^ 100 GeV and confirm the two- 



Ultimate synchrotron cutoff 5 

componentness in this way even without detection of the TeV emis- 
sion by large ground-based detectors. 

Although the possibility of the extreme acceleration up to the 
energy 7ni in internal shocks is not evident, it is not yet ruled out. 
Do independent signatures exist of the photon breding mechanism? 
A possible one is the correlation between the disc brightness and 
the spectral shape like those in Fig. (2] 

If the total power of the jet is known, then the high radiative 
efficiency, i.e. the ratio of the gamma-ray luminosity to the total jet 
power, of e.g. more than 10 per cent can give an additional sup- 
port to the photon breeding hypothesis and argue against the inter- 
nal shock scenario. In the internal shock mechanism the luminosity 
budget is limited by the internal energy of the jet, while in the case 
of the photon breeding a large fraction of the total bulk energy of 
the jet can be converted into radiation. 

Another possible clue is the compactness threshold: at a low 
compactness of the disc radiation, the photon breeding does not 
work (the threshold for model C is L^/Ris ~ (1 ~ 10) x 
10*^ erg s~^). Usually it is difficult to estimate the distance R from 
the disc to the emission region, however, sometimes approximate 
constraints can be obtained using the time variability or the mass 
of the black hole (R can not be smaller than a few i?g). In this 
way the photon breeding mechanism can be revealed statistically: 
sources with the low estimated compactness should be less effi- 
cient and have different spectra. When the photon breeding turns 
on, the luminosity should increase sharply. Therefore one can ex- 
pect a kind of bimodality in the efficiency and probably in the high- 
energy spectra of blazars. 
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